Journal  of  Power  Sources  234  (2013)  69-81 


ELSEVIER 


Contents  lists  available  at  SciVerse  ScienceDirect 

Journal  of  Power  Sources 

journal  homepage:  www.elsevier.com/locate/jpowsour 


Charge  transfer  and  formation  of  Ce3+  upon  adsorption  of  metal  atom 
M  (M  =  Cu,  Ag,  Au)  on  Ce02  (100)  surface 

Li-Jiang  Chen3,  Yuanhao Tang3,  Lixia  Cui3,  Chuying  0uyangc,  Siqi  Shi3,b’* 

1  School  of  Sciences,  Zhejiang  Sci-Tech  University,  Xiasha  College  Park,  Hangzhou  310018,  China 
b  School  of  Materials  Science  and  Engineering,  Shanghai  University,  Shanghai  200444,  China 
c  Department  of  Physics,  Jiangxi  Normal  University,  Nanchang  330022,  China 


HIGHLIGHTS 


►  Adsorption  of  M  (M  =  Cu,  Ag,  Au)  on  Ce02(100)  causes  in  charge  transfer  from  M  to  Ce. 

►  For  the  same  adsorption  site,  there  exist  multiple  adsorption  configurations. 

►  Our  theoretical  results  are  compared  wherever  possible  with  experiment. 
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Adsorption  behavior  of  metal  atom  M  [M  =  Cu,  Ag,  Au)  on  CeO2(100),  a  technologically  important  cat¬ 
alytic  support  surface,  is  studied  by  using  first-principles  density-functional  theory  calculations  with  the 
on-site  Coulomb  interaction  taken  into  account  (DFT  +  U).  Prior  to  geometry  optimization,  small  dis¬ 
tortions  in  selected  Ce-0  distances  are  imposed  to  explore  the  energetics  associated  with  reduction  of 
Ce4+  to  Ce3+  due  to  charge  transfer  to  Ce  during  M  adsorption.  Results  indicate  that  the  adsorption  is 
accompanied  by  an  electron  charge  transfer  between  neutral  metal  atom  and  neighboring  Ce4+  cation. 
For  the  same  adsorption  site,  there  exist  multiple  adsorption  configurations  with  the  excess  electron 
localized  at  various  nearest  neighboring  Ce  ion  relative  to  the  M.  The  order  on  adsorption  energies  of  M 
on  CeO2(100)  surface  follows:  Cu  >  Au  >  Ag.  The  atomic  and  electronic  structures  of  M  on  Ce02(100)  are 
compared  with  those  of  M  on  CeC>2(lll)  and  CeO2(110). 

©  2013  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Materials  based  on  CeC>2  have  attracted  much  attention  because 
of  their  uses  in  a  variety  of  processes  such  as  automotive  exhaust 
catalysis,  water-gas  shift  (WGS)  reactions,  fuel  cells  and  production 
or  purification  of  hydrogen  [1,2].  These  applications  of  CeC>2  are 
largely  attributed  to  its  remarkable  oxygen-storage  capability  and 
the  strong  interaction  between  loaded  metal  atom  and  Ce02  sup¬ 
port  [3].  This  has  been  confirmed  by  examining  the  oxidation  of  CO 
[4],  reduction  of  N0X  [5]  and  destruction  of  S02  [6]  on  metal 
nanoparticles  dispersed  on  Ce02.  Some  investigations  [7—13]  have 
shown  that  the  dispersion  of  M  (M  =  Cu,  Au,  Pt,  Rh)  in  a  Ce02  matrix 
can  produce  materials  with  novel  structural  and/or  electronic 
properties,  especially  the  formation  of  the  ionic  M—O  bonds, 
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playing  an  important  role  in  the  conversion  of  CO  (NO  or  SO2)  in 
M/Ce O2  catalysts.  In  other  words,  the  kinetic  behaviors  of  the  Ce02 
supported  catalysts  depend  strongly  on  the  morphology  and 
chemical  state  of  the  active  metal  species.  For  example,  Kydd  et  al. 
[7]  and  Avgouropoulos  et  al.  [13]  have  observed  the  conversion  of 
CO  in  the  Cu/Ce02  system  with  Cu— O  or  Cu— 0— Ce  bonds  by  Flame 
Spray  Pyrolysis  (FSP),  Brunaner— Emmett— Teller  (BET)  adsorption 
isotherms,  X-ray  diffraction  (XRD),  High  Resolution  Transmission 
Electron  Microscopy  (HRTEM),  Temperature  Programmed  Reduc¬ 
tion  (TPR)  and  Temperature  Programmed  Desorption  (TPD). 

Performance  of  Ce02  supported  catalysts  is  also  relevant  to  the 
electronic  structures  of  the  adsorbed  metal  atom.  It  is  generally 
recognized  that  the  adsorbed  metal  atom  is  oxidized  and  thus  brings 
in  transformation  of  Ce4+  — >  Ce3+.  The  adsorbed  metal  ions  are  active 
in  the  oxidation  catalytic  reaction.  Some  theoretical  investigations 
have  revealed  that  the  metal  atoms  could  form  relatively  strong 
bonds  with  oxygen  atoms  by  transferring  their  outer  electron  to  the 
substrate  and  thus  Ce02  surface  is  reduced.  For  example,  Branda  et  al. 
[14]  andTang  etal.  [15]  have  investigated  the  adsorption  behavior  of 
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Cu,  Ag  and  Au  on  Ce02(lll)  surface  and  found  JW+(M  =  Cu,  Ag,  Au). 
Hernandez  etal.  [16]  suggested  that  the  high  activity  of  Au/Ce02(lll ) 
system  should  be  due  to  Au  ions  resulting  from  charge  transfer  from 
Au  to  Ce  atom.  Therefore,  it  is  of  great  importance  to  clarify  the  effect 
of  the  adsorption  of  metal  atom  on  its  electronic  states  or  the  charge 
redistribution  at  the  M/Ce O2  interfaces. 

Studies  on  atomic  and  electronic  structure  of  metal  atom  M 
(M  =  Cu,  Ag,  Au)  on  the  Ce02(lll)  and  (110)  surfaces  can  be  found 
in  Refs.  [14—20].  In  our  preceding  papers  [15,17],  we  have  made 
detailed  analyses  on  the  adsorption  behaviors  of  metal  M  (M  =  Cu, 
Ag,  Au)  atom  on  CeC>2(lll)  and  (110)  surfaces  by  using  first- 
principles  calculations  based  on  density  functional  theory  (DFT). 
Charge  redistributions  with  electron  localized  on  various  Ce  cations 
neighboring  to  adsorbed  atom  M  (M  =  Cu,  Ag,  Au)  are  observed  and 
can  be  attributed  to  local  structural  distortion  effects.  We  demon¬ 
strate  that  the  electron  localized  on  the  Ce  ion  away  from  the  metal 
M  (M  =  Cu,  Ag,  Au)  atom  and  the  adatom  is  oxidized  to  be  M+  on 
CeC^lll)  and  CeC^llO)  surfaces,  and  Cu  adatom  on  CeC>2(110) 
surface  can  also  be  oxidized  to  Cu2+.  More  importantly,  previous 
controversies  about  the  valence  state  of  the  adsorbed  Au  atom  on 
Ce02(lll )  surface  have  been  successfully  clarified  via  analysis  of  the 
local  structural  distortions  [15].  Nevertheless,  a  careful  examination 
of  the  literature  revealed  no  other  attempts  to  use  first-principles 
calculation  to  study  the  adsorption  behavior  of  M  adsorption  on 
CeC>2(100).  Following  the  classification  by  Sayle  et  al.  [21], 
Ce02(100)  surface  can  be  described  as  the  kind  of  one,  on  which 
exists  a  net  dipole  normal  to  the  surface.  Such  a  surface  would  be 
described  as  a  polar  one  and  predicted  to  have  a  high  surface  en¬ 
ergy  or  catalytic  activity  [22—31  ].  Besides,  several  studies  [24,31,32] 
have  indicated  that  the  (100)  surface  undergoes  a  drastic  structure 
rearrangement  by  usage  of  scanning  tunneling  microscope  (STM), 
low  energy  electron  diffraction  (LEED)  and  the  molecular  dynamics 
simulation.  Accordingly,  compared  with  the  non-polar  (111)  and 
(110)  surfaces,  the  polar  (100)  surface  seems  more  reactive  and 
consequently  more  important  to  comprehensively  understand  the 
catalytic  behaviors  of  M/CeC>2  (M  =  Cu,  Ag,  Au)  systems.  For 


example,  experimental  works  of  Zhou  et  al.  [33]  and  Aneggi  et  al. 
[34]  have  provided  an  evidence  of  the  surface  sensitivity  of  the  CO— 
Ce02  interaction  by  studying  the  catalytic  activity  of  Ce02  nanorods 
and  nanoparticles  and  polycrystalline  Ce02;  it  was  demonstrated 
that  Ce02  nanorods,  with  well-defined  reactive  planes  (001)  and 

(110) ,  are  more  active  than  “classical”  ceria  nanoparticles  with  the 

(111 )  planes.  Trtik  et  al.  [35]  prepared  Ce02(001 )  epitaxial  thin  films 
deposited  on  Si(001 )  with  yttria-stabilized  zirconia  buffers  at  room 
temperature  by  pulsed-laser  deposition.  Therefore,  it  is  worthwhile 
to  investigate  the  relative  stabilities  and  electronic  redistribution  of 
the  adsorption  of  Cu,  Ag  or  Au  on  the  (100)  surface  of  ceria  and  to 
ascertain  whether  the  polarity  and  instability  of  the  CeO2(100) 
surface  would  in  turn  affect  the  adsorption  structures  or  electronic 
properties  of  metal  atom  on  CeC>2  surfaces. 

In  this  paper,  we  examine  the  effect  of  the  adsorption  of  M 
(M  =  Cu,  Ag,  Au)  atom  on  the  Ce-4/electron  localization  on  CeO2(100) 
surface  and  the  associated  structural  stability  by  using  the  first- 
principles  calculations  based  on  the  density  functional  theory 
(DFT).  We  find  the  most  energetically  favorable  structure  is  M 
adsorption  on  O-bridge  (Ou— Ou)  site.  Charge  redistributions,  with 
electrons  localized  on  INN,  2NN  and  3NN  (The  INN,  2NN  and  3NN 
notations  represent  the  Ce  ions  which  are  the  first,  second  and  third 
nearest  neighbors  to  the  adatoms,  respectively.)  Ce  ions  to  M  are 
observed  with  Ce4+  reduced  to  Ce3+  in  most  cases.  In  addition, 
a  detailed  comparison  of  the  atomic  and  electronic  structure  of  M  on 
Ce02(lll )  and  (110)  surfaces  is  made  for  the  selected  adsorption  sites. 
An  important  difference  from  our  previous  studies  on  M  (M  =  Cu,  Ag, 
Au)  on  CeC>2(  111 )  and  (110)  surfaces  [  15,17  ]  is  the  observation  that  the 
M  adsorption  configuration  with  the  excess  electron  located  at  the 
first  nearest  neighboring  Ce4+  ion  is  more  energetically  favorable.  The 
present  methodology  and  results  provide  a  useful  basis  of  compari¬ 
son  for  experimental  studies  focused  on  catalyst  designs. 

The  computational  details  are  outlined  in  Section  2.  And  Section 
3  details  the  adsorption  models,  energetics,  and  discussions  of  the 
major  theoretical  predictions.  The  most  significant  conclusions 
from  this  study  are  listed  in  Section  4. 
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Fig.  1.  CeO2(100)  p( 2  x  2)  surface  and  M  (M  =  Cu,  Ag,  Au)  adsorption  sites  (labeled  as  “x”).  (a)  Oblique  view  along  [001];  (b)  View  of  surface  along  [100].  Red  and  blue  spheres 
represent  O  and  Ce,  respectively.  Large  spheres  denote  the  first  layer  ions.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version 
of  this  article.) 
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2.  Computational  details  and  supercell  construction 

The  fully  spin-polarized  calculations  are  performed  by  using  the 
frozen-core  all-electron  projected  augmented  wave  (PAW)  method 

[36] ,  as  implemented  in  the  ab  initio  total  energy  and  molecular 
dynamic  program  VASP  (the  Vienna  ab  inito  simulation  package) 

[37] .  The  Perdew— Burke— Ernzerhof  (PBE)  functional  [38]  is  used 
for  the  exchange  correlation.  To  account  for  strong  on-site  Coulomb 
repulsion  among  the  Ce  4/  electrons,  a  Hubbard  parameter  U  is 
added  to  the  PBE  functional  with  the  rotationally  invariant 
approach  formulated  by  Dudarev  et  al.  [39],  in  which  only  the 
difference  (t/eff  =  U  -  J)  between  the  Coulomb  U  and  exchange  J 
parameters  must  be  specified.  As  suggested  by  Fabris  et  al.  [40]  for 


CeC>2  using  the  linear-response  approach  of  Cococcioni  and  de 
Gironcoli  [41],  the  Ueff  value  of  4.5  eV  is  used.  The  choice  of  this 
methodology  (including  PBE  functional,  rotationally  invariant 
approach  and  Ueff  value)  is  supported  by  previous  studies  on  CeC>2 
systems  that  have  shown  a  very  good  agreement  with  available 
results  for  bulk  CeC>2  and  Ce02(lll)  and  (110)  surfaces  [15,17,42— 
45],  The  present  PBE  +  U  calculations  give  an  equilibrium  lattice 
parameter  of  5.49  A  [46],  which  is  close  to  the  experimental  value  of 
5.41  A  [47,48], 

Valence  electron  configurations  for  the  elemental  constituents 
are  as  follows:  Ce-5s25p66s25d14/1,  0-2s22p4,  Cu-Bd’^s1,  Ag- 
4d105s1  and  Au-Sd^s1.  The  Monkhorst— Pack  scheme  [49]  based 
on  the  2x2x1  k-point  grid  and  the  Gaussian  smearing  of  0.20  eV 


(c)  Ce-top 


Fig.  2.  Structural  models  of  M  (M  =  Cu,  Ag,  Au)  atoms  adsorbed  on  Ce02(100).  (a)  Ou  (O-top),  (b)  Ou-Ou  (O-bridge),  (c)  Ce-top.  Red,  blue  and  yellow  spheres  denote  0,  Ce  and  M 
atoms,  respectively.  The  INN  to  3NN  notations  indicate  the  Ce  ions  which  are  first,  second  and  third  nearest  neighbors,  respectively,  to  the  adsorbed  M  atom.  (For  interpretation  of 
the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 
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Table  1 

Adsorption  energies,  bond  lengths  and  magnetic  spin  moments  of  metal  atoms  for  the  adsorption  configurations  of  Cu  atom  on  CeO2(100).  For  a  given  adsorption  site, 
configuration  numbers  (e.g.  conf.  2,  INN)  are  associated  with  the  distortion  applied  to  Ce— 0  distances  listed  as  dCe- o.  which  are  the  bond  lengths  between  the  Ce3+  and  the 
nearest  neighboring  0  ions.  The  undistorted  dCe- o  is  2.38  A.  Eads  is  the  computed  adsorption  energy,  dCu-o  is  the  distance  between  the  adsorbed  Cu  and  the  nearest  neighboring 
0  ions  (an  0  ion  for  the  Ou  site  and  two  0  ions  for  the  Ou— Ou  and  Ce-top  sites),  and  dcu-ce  are  distances  from  the  adsorbed  Cu  to  INN,  2NN,  and  3NN  Ce  ions  (INN,  the  first 
nearest  neighbor,  2NN,  the  second  nearest  neighbor  and  3NN,  the  third  nearest  neighbor  Ce  ion  from  the  Cu  adatom),  respectively,  following  geometry  optimization.  Note  that 
integers  in  dCu-o  (A)  and  dCu-ce  (A)  denote  the  number  of  0  and  Ce  ions,  respectively,  within  the  distances  relative  to  the  adsorbed  Cu  ion.  Magnetic  spin  moments  of  reduced 
Ce  (in  bold)  and  oxidized  Cu  ions  are,  respectively,  fiCe  and  iiCu- 


Site 

Eads(eV) 

dcu-o  (A) 

dc u— Ce  (A) 

PCeiPB) 

dee- O  (A) 

MCu(Mb) 

Ou  (conf.  1) 

2.093 

1  x  1.772 

1  x  3.667 

0.526 

1  X  2.429 

0.000 

1  x  3.667 

0.527 

2  x  2.272 

4  x  5.298 

0.000 

1  x  2.566 

2  x  6.088 

0.000 

1  x  2.236 

2  x  6.603 

0.000 

1  x  2.422 

Ou  (conf.  2,  INN) 

2.376 

1  x  1.773 

1  x  3.584 

0.000 

1  x  2.478 

0.000 

1  x  3.783 

0.958 

1  x  2.326 

1  x  5.216 

0.000 

1  x  2.321 

1  x  5.249 

0.000 

1  x  2.662 

1  x  5.385 

0.000 

1  x  2.358 

1  x  5.421 

0.000 

1  x  2.529 

1  x  6.111 

0.000 

1  x  6.125 

0.000 

1  x  6.363 

0.000 

1  x  6.855 

0.000 

Ou  (conf.  3,  2NN) 

2.296 

1  x  1.768 

1  x  3.630 

0.000 

1  x  2.589 

0.000 

1  x  3.656 

0.000 

1  x  2.332 

1  x  5.204 

0.000 

1  x  2.345 

1  x  5.215 

0.000 

1  x  2.592 

1  x  5.245 

0.000 

1  x  2.423 

1  x  5.176 

0.957 

1  x  2.424 

1  x  6.009 

0.000 

1  x  6.000 

0.000 

1  x  6.560 

0.000 

1  x  6.537 

0.000 

Ou  (conf.  4,  3NN) 

2.316 

1  x  1.772 

1  x  3.611 

0.000 

1  x  2.604 

0.000 

1  x  3.663 

0.000 

1  x  2.332 

2  x  5.253 

0.000 

1  x  2.337 

1  x  5.282 

0.000 

1  x  2.550 

1  x  5.285 

0.000 

1  x  2.442 

1  x  6.063 

0.000 

1  x  2.396 

1  x  6.064 

0.000 

1  x  6.522 

0.963 

1  x  6.549 

0.000 

Ou-Ou  (conf.  1 ) 

3.702 

2  x  1.810 

2  x  3.182 

2  x  0.288 

1  x  2.452 

0.000 

2  x  3.182 

2  x  0.287 

1  x  2.251 

1  x  4.249 

0.000 

1  x  2.262 

2  x  5.744 

0.000 

1  x  2.539 

2  x  5.786 

0.000 

1  x  2.197 

4  x  6.312 

0.000 

1  x  2.378 

Ou-Ou  (conf.  2,  INN) 

4.079 

1  x  1.809 

1  x  3.149 

0.000 

1  x  2.491 

0.000 

1  x  1.813 

1  x  3.158 

0.000 

1  x  2.319 

1  x  3.225 

0.000 

1  x  2.347 

1  x  3.254 

0.958 

1  x  2.674 

1  x  4.259 

0.000 

1  x  2.349 

1  x  5.748 

0.000 

1  x  2.572 

2  x  5.762 

0.000 

1  x  5.843 

0.000 

1  x  6.242 

0.000 

1  x  6.249 

0.000 

1  x  6.286 

0.000 

1  x  6.311 

0.000 

Ou-Ou  (conf.  3,  2NN) 

3.916 

1  x  1.809 

1  x  3.193 

0.000 

1  x  2.484 

0.000 

1  x  1.810 

1  x  3.197 

0.000 

1  x  2.421 

1  x  3.199 

0.000 

1  x  2.419 

1  x  3.207 

0.000 

1  x  2.486 

1  x  4.264 

0.859 

1  x  2.403 

1  x  5.776 

0.000 

1  x  2.576 

1  x  5.786 

0.000 

1  x  2.586 

1  x  5.811 

0.000 

1  x  2.403 

1  x  5.812 

0.000 

1  x  6.332 

0.000 

1  x  6.316 

0.000 

1  x  6.328 

0.000 

1  x  6.319 

0.000 

Ou-Ou  (conf.  4,  3NN) 

3.920 

1  x  1.815 

1  x  3.178 

0.000 

1  x  2.573 

0.000 

1  x  1.813 

1  x  3.184 

0.000 

1  x  2.332 

1  x  3.186 

0.000 

1  x  2.333 
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Table  1  (continued) 


Site 

Fads(eV) 

d cu-o  (A) 

dcu-Ce  (A) 

MCe(MB) 

dee- 

-o(A) 

MCu(Mb) 

1 

X  3.207 

0.000 

1  X 

2.643 

1 

x  4.261 

0.000 

1  X 

2.378 

1 

x  5.743 

0.000 

1  X 

2.392 

1 

x  5.783 

0.000 

1 

x  5.790 

0.000 

1 

x  5.800 

0.000 

1 

x  6.244 

0.909 

1 

x  6.326 

0.000 

1 

x  6.266 

0.000 

1 

x  6.347 

0.000 

Ce-top  (conf.  1) 

2.221 

1  x  1.948 

1 

x  2.444 

0.000 

1  X 

2.530 

0.000 

1  x  1.948 

1 

x  4.510 

0.000 

1  X 

2.237 

1 

x  4.599 

0.000 

1  X 

2.243 

1 

x  4.602 

0.600 

1  X 

2.469 

1 

x  4.635 

0.600 

1  X 

2.220 

1  X 

2.192 

Ce-top  (conf.  2,  2NN) 

2.638 

1  x  1.897 

1 

x  2.710 

0.000 

1  X 

2.633 

0.000 

1  x  1.907 

1 

x  4.078 

0.000 

1  X 

2.304 

1 

x  4.632 

0.000 

1  X 

2.311 

1 

x  4.675 

0.000 

1  X 

2.378 

1 

x  5.370 

0.969 

1  X 

2.293 

1  X 

3.463 

Distances  between  Cu  and  reduced  Ce  are  labelled  in  bold. 

are  employed  for  the  Brillouin-Zone  integrals.  The  cutoff  energy  for 

sites  must  be  considered.  (For  interpretation  of  the  references  to 

the  plane-wave  is  400  eV.  This  set  of  parameters  assures  a  total 
energy  convergence  of  10  7  eV  unit  cell-1.  The  structure  optimi¬ 
zations  are  performed  until  all  the  Hellmann— Feynman  forces  are 
converged  to  be  less  than  0.01  eV  A-1.  For  a  more  accurate  calcu¬ 
lation  of  the  electronic  density  of  states  (DOS),  we  apply  4x4x11;- 
point  grid  and  the  modified  tetrahedron  method  [50]. 

The  (100)  surface  is  modeled  as  an  oxygen-terminated  p( 2  x  2) 
( p  =  primitive)  slab  (red  enlarged  spheres  denote  surface  O  ions  in 
Fig.  1(a))  with  half  of  the  surface  oxygen  atoms  moved  from  one 
side  of  the  slab  to  the  other  in  a  checkerboard  style  to  fulfill  the 
Ce02  formula  and  to  avoid  a  dipole  moment  normal  to  the  surface 
[32].  Fig.  1(a)  shows  the  slab  has  a  26.98  A  thickness  and  contains 
32  CeC>2  formula  units  (nine  atomic  layers)  and  a  vacuum  layer  of 
16  A.  The  x  [010]  andy  [001  ]  dimensions  of  the  slab  are  both  fixed  at 
the  10.98  A  as  shown  in  Fig.  1(b),  which  is  based  on  the  optimized 
5.49  A  Ce02  lattice  parameter.  Although  the  relaxation  in  the 
CeO2(100)  surface  is  complex  (because  of  the  removal  of  half 
a  monolayer  of  outer  oxygen  atoms)  as  compared  with  the  (111) 
and  (110)  surfaces,  surface  energy  convergence  tests  [23]  show  that 
relaxation  in  the  Ce02(100)  surface  occurs  within  the  initial  five 


color  in  this  paragraph,  the  reader  is  referred  to  the  web  version  of 
this  article.) 

Prior  to  geometry  optimization,  small  structural  distortions  were 
introduced  into  the  as-built  models  to  achieve  various  localized 
solutions.  This  required  slight  adjustment  of  distances  between 
a  Ce4+  ion  and  its  nearest-neighboring  O  ions  to  be  consistent  with 
those  of  a  Ce3+  ion.  The  positions  of  M  on  the  surface  and  the 
associated  bond  length  adjustments  resulted  in  a  range  of  model 
configurations  which  we  now  explore.  Experimentally,  these  slight 
distortions  can  be  realized  by  introducing  impurities  or  manipu¬ 
lating  the  crystal  morphology  during  catalytic  material  preparation 
[51,52].  The  adsorption  energy,  £ads.  of  M  atom  on  Ce02(100)  is 
calculated  via: 


£ads  =  £ce02  +  “  ^M/CeOj 


(1) 


where  £Ceo2.  Em  and  £M/ce02  correspond  to  the  calculated  total 
energies  of  the  clean  Ce02(100)p(2  x  2),  isolated  M  atom  and 
combined  M/Ce02(100)p(2  x  2)  system,  respectively. 


layers  and  the  atomic  structures  at  the  center  of  the  slabs  are  bulk¬ 
like.  Our  calculated  results  show  that  the  interlayer  distance  in  the 
Ce  sublattice  increases  slightly  (expanded  ca.  0.04  A  near  the  sur¬ 
face).  The  original  coplanar  O  atoms  in  the  third  and  fifth  atomic 
layers  are  split  into  two  sets  with  the  separation  of  0.46  and  0.17  A 
respectively.  The  bottom  four  layers  are  fixed  at  their  bulk  positions 
(see  region  enclosed  by  dashed  line  in  Fig.  1(a))  and  remaining 
layers  are  allowed  to  relax.  (For  interpretation  of  the  references  to 
color  in  this  paragraph,  the  reader  is  referred  to  the  web  version  of 
this  article.) 

Fig.  1(b)  shows  three  adsorption  positions  of  M  (M  =  Cu,  Ag,  Au) 
atoms  on  the  Ce02(100)  surface.  Here,  the  following  definitions 
apply:  Ou  (a  surface  O-top  site);  Ou— Ou  (a  surface  O-bridge  site), 
and  Ce-top  (a  surface  Ce-top  site).  Each  occurrence  of  x  in  Fig.  1(b) 
denotes  an  M  adsorption  site  on  CeO2(100).  Fig.  2  shows  the  posi¬ 
tion  of  an  adsorbed  M  (yellow)  at  Ou,  directly  above  a  surface  O  ion 
(Fig.  2(a)),  0LI— Ou,  an  oxygen  ion  bridge  site  (Fig.  2(b)),  and  Ce-top, 
a  Ce  ion  top  site  (Fig.  2(c))  on  CeO2(100)  surface.  This  is  different 
from  previous  studies  ofM/Ce02(lll)  and  M/Ce O2(110),  where  four 
(Od,  Ou,  Ou-Ou,  Ce— Ce)  and  three  (Ou,  Ou— Ou,  Ce— Ce)  adsorption 


3.  Results  and  discussion 

Identified  by  computing  £ads  ate  the  Ce3+  location,  the  M— O,  M— 
Ce  and  Ce— O  bond  lengths,  magnetic  spin  moments  of  M  and  Ce, 
and  various  configurations  for  4/electrons  localized  at  neighboring 
Ce  ions  to  the  M  (M  =  Cu,  Ag,  Au).  These  quantities  are  listed  in 
Tables  1—3  for  Cu,  Ag  and  Au,  respectively.  Results  without  small 
structural  distortions  at  three  sites  are  also  listed,  although  the 
following  discussion  only  refers  to  the  distorted  cases  unless  oth¬ 
erwise  explicitly  stated.  We  consider  the  Ou— Ou  adsorption  site 
(shown  in  Fig.  2(b))  as  an  example.  There  are  three  different  dis¬ 
tances  of  Ce3+  ions  around  the  adsorbed  M  ion.  These  are  denoted 
as  INN,  2NN  and  3NN,  respectively,  as  mentioned  above.  The 
computed  spin  moments  of  Ce  and  Cu  in  Table  1  show  that  the  s 
electron  from  the  Cu  atom  can  localize  on  neighboring  Ce  ions  since 
Cu-3d10  has  no  moment  associated  with  it.  Note  that  for  the  Ou— Ou 
site,  the  2NN  Ce  ion  is  located  in  the  subsurface  (i.e.,  the  fourth 
atomic  layer  in  Fig.  1(a))  whereas  the  INN  and  3NN  Ce  ions  are  at 
the  surface  (i.e.,  the  second  atomic  layer  in  Fig.  1(a)).  This  is 
opposite  to  the  case  of  the  Ou— Ou  site  for  JVf/CeO2(110),  and  is 
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Table  2 

Adsorption  energies,  bond  lengths  and  magnetic  spin  moments  of  metal  atoms  for  the  adsorption  configurations  of  Ag  atom  on  CeO2(100).  For  a  given  adsorption  site, 
configuration  numbers  (e.g.  conf.  2,  INN)  are  associated  with  the  distortion  applied  to  Ce— 0  distances  listed  as  dCe- o.  which  are  the  bond  lengths  between  the  Ce3+  and  the 
nearest  neighboring  0  ions.  The  undistorted  dCe- o  is  2.38  A.  Eads  is  the  computed  adsorption  energy,  dAg_0  is  the  distance  between  the  adsorbed  Ag  and  the  nearest  neighboring 


0  ions  (an  0  ion  for  the  Ou  site  and  two  0  ions  for  the  Ou— Ou  and  Ce-top  sites),  and  dAg_ce  are  distances  from  the  adsorbed  Ag  to  INN,  2NN,  and  3NN  Ce  ions  (INN,  the  first 
nearest  neighbor,  2NN,  the  second  nearest  neighbor  and  3NN,  the  third  nearest  neighbor  Ce  ion  from  the  Ag  adatom),  respectively,  following  geometry  optimization.  Note  that 
integers  in  dAg_G  (A)  and  dAg_Ce  (A)  denote  the  number  of  0  and  Ce  ions,  respectively,  within  the  distances  relative  to  the  adsorbed  Ag  ion.  Magnetic  spin  moments  of  reduced 
Ce  (in  bold)  and  oxidized  Ag  ions  are,  respectively,  nCe  and  ^Ag. 

Site  E^s  (eV)  dAg-o  (A) 

dAg-Ce  (A) 

PCeiPB) 

dce-0  (A) 

MAg(^B) 

Ou  (conf.  1)  1.307  1  x  2.059 

1  x  3.858 

0.496 

1  X  2.448 

0.000 

1  x  3.858 

0.495 

2  X  2.272 

4  x  5.440 

0.000 

1  x  2.557 

2  x  6.324 

0.000 

1  x  2.247 

2  x  6.711 

0.000 

1  x  2.379 

Ou  (conf.  2,  INN)  1.437  1  x  2.062 

1  x  3.703 

0.000 

1  x  2.486 

0.000 

1  x  4.057 

0.901 

2  x  2.307 

1  x  5.315 

0.000 

1  x  2.642 

1  x  5.310 

0.000 

1  x  2.356 

1  x  5.595 

0.000 

1  x  2.506 

1  x  5.589 

0.000 

1  x  6.341 

0.000 

1  x  6.338 

0.000 

1  x  6.320 

0.000 

1  x  7.108 

0.000 

Ou  (conf.  3,  2NN)  1.369  1  x  2.067 

1  x  3.879 

0.000 

1  x  2.584 

0.000 

1  x  3.862 

0.000 

1  x  2.323 

1  x  5.360 

0.887 

1  x  2.326 

1  x  5.371 

0.000 

1  x  2.589 

1  x  5.406 

0.000 

1  x  2.402 

1  x  5.402 

0.000 

1  x  2.401 

1  x  6.299 

0.000 

1  x  6.291 

0.000 

1  x  6.708 

0.000 

1  x  6.699 

0.000 

Ou  (conf.  4,  3NN)  1.561  1  x  2.066 

1  x  3.842 

0.000 

1  x  2.592 

0.000 

1  x  3.880 

0.000 

1  x  2.333 

1  x  5.447 

0.000 

1  x  2.328 

1  x  5.446 

0.000 

1  x  2.567 

1  x  5.423 

0.000 

1  x  2.446 

1  x  5.420 

0.000 

1  x  2.41 1 

1  x  6.641 

0.962 

2  x  6.336 

0.000 

1  x  6.666 

0.000 

Ou-Ou  (conf.  1)  2.045  2  x  2.077 

2  x  3.325 

2  x  0.263 

1  x  2.450 

0.000 

2  x  3.325 

2  x  0.264 

1  x  2.245 

1  x  4.501 

0.000 

1  x  2.264 

2  x  5.950 

0.000 

1  x  2.534 

2  x  5.975 

0.000 

1  x  2.204 

4  x  6.383 

0.000 

1  x  2.358 

Ou-Ou  (conf.  2,  INN)  2.437  2  x  2.070 

1  x  3.301 

0.000 

1  x  2.499 

0.000 

1  x  3.310 

0.000 

1  x  2.323 

1  x  3.378 

0.000 

1  x  2.339 

1  x  3.378 

0.960 

1  x  2.638 

1  x  4.532 

0.000 

1  x  2.362 

1  x  5.952 

0.000 

1  x  2.557 

1  x  5.970 

0.000 

1  x  5.984 

0.000 

1  x  6.041 

0.000 

1  x  6.325 

0.000 

1  x  6.404 

0.000 

1  x  6.354 

0.000 

1  x  6.321 

0.000 

Ou-0,,  (conf.  3,  2NN)  2.263  2  x  2.073 

1  x  3.352 

0.000 

1  x  2.490 

0.000 

1  x  3.354 

0.000 

1  x  2.417 

1  x  3.341 

0.000 

1  x  2.418 

1  x  3.339 

0.000 

1  x  2.492 

1  x  4.527 

0.848 

1  x  2.394 

2  x  6.011 

0.000 

1  x  2.577 

1  x  5.991 

0.000 

1  x  2.578 

1  x  5.993 

0.000 

1  x  2.396 

Ou-Ou  (conf.  4.  3NN)  2.205  1  x  2.077 

1  x  3.351 

0.000 

1  x  2.557 

0.000 

1  x  2.078 

1  x  3.338 

0.000 

1  x  2.327 

1  x  3.347 

0.000 

1  x  2.324 

1  x  3.324 

0.000 

1  x  2.625 

1  x  4.539 

0.000 

1  x  2.368 

1  x  5.981 

0.000 

1  x  2.391 

1  x  5.965 

0.000 
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Table  2  (continued) 


Site 

Fads  (eV) 

dAg-o  (A) 

dAg-Ce  (A) 

MCe(^B) 

dce-o  (A) 

/Wmb) 

1  x  6.009 

0.000 

1  x  6.003 

0.000 

1  x  6.422 

0.000 

1  x  6.331 

0.000 

1  x  6.324 

0.861 

1  x  6.395 

0.000 

Ce-top  (conf.  1 ) 

1.358 

2  x  2.316 

1  x  2.927 

0.000 

1  x  2.548 

0.000 

2  x  4.749 

0.000 

2  x  2.276 

2  x  4.730 

2  x  0.501 

1  x  2.541 

1  x  2.452 

1  x  2.225 

Ce-top  (conf.  2,  2NN) 

1.670 

1  x  2.293 

1  x  2.944 

0.000 

1  x  2.534 

0.000 

1  x  2.326 

1  x  4.912 

0.972 

2  x  2.321 

2  x  4.774 

0.000 

1  x  2.595 

1  x  4.599 

0.000 

1  x  2.387 

1  x  2.661 

Distances  between  Ag  and  reduced  Ce  are  labelled  in  bold. 

different  from  the  case  of  the  Ou— Ou  site  for  M/Ce02(in ),  in  which 
the  INN,  2NN  and  3NN  Ce  ions  are  at  the  surface.  All  three  con¬ 
sidered  sites  would  provide  different  configurations  with  the 
excess  electron  localized  on  Ce  ions  which  are  different  neighbors 
to  the  adsorbed  atoms.  For  the  convenience  of  comparison,  the 
delocalized  configurations,  which  are  less  stable  than  the  localized 
ones,  are  also  listed  in  Tables  1—3,  although  the  configuration 
reported  in  this  work  refers  to  the  localized  one  unless  otherwise 
explicitly  stated. 

3.1.  Cu/CeO2(100) 

Table  1  summarizes  the  adsorption  energies,  £ads  (eV),  from 
Eq.  (1),  relevant  bond  lengths  and  magnetic  spin  moments  for  all 
the  possible  adsorption  configurations  of  Cu  adsorbed  on  CeC>2 
(100)  surface,  dcu-ce  is  the  distance  (A)  from  adsorbed  Cu  to  INN, 
2NN  or  3NN  Ce  ion;  dcu-o  is  the  distance  between  the  adsorbed  Cu 
and  the  nearest  neighboring  O  ion;  dce-o  (A)  is  the  Ce3+— O  bond 
length;  jice  and  ^cu  are  spin  moments  (jlb)  of  Ce  and  Cu,  respec¬ 
tively.  Integers  that  precede  dcu-o  and  dcu-Ce  denote  the  numbers 
of  NN  0  and  Ce  ions  relative  to  the  adsorbed  Cu  atom,  respectively. 
We  can  clearly  see  that,  among  the  three  types  of  adsorption  sites, 
the  adsorption  of  Cu  atom  at  the  two-fold  Ou— Ou  site  exhibits  the 
largest  adsorption  energies  (3.916,  3.920  and  4.079  eV  depending 
on  the  location  of  Ce3+  relative  to  the  Cu  atom).  The  second  stable 
one  is  the  adsorption  at  the  Ce-top  site  with  two-fold  coordinated 
Cu— O  bonds  (2.638  eV).  The  adsorption  energies  (2.296,  2.316  and 
2.376  eV)  at  the  Ou  site  are  slightly  smaller  than  that  at  the  Ce-top 
site.  On  the  other  hand.  Note  that  difference  between  adsorption 
energies  of  any  two  configurations  at  the  Ou  or  Ou— Ou  site  is  rela¬ 
tively  minor  (0.004—0.163  eV)  and  that  a  distribution  of  these 
configurations  could  be  present  in  the  experiments.  Distribution  of 
Ti3+  (Ce3+)  sites  on  the  reduced  TiO2(110)  (Ce02(lll))  surface  is 
a  good  example  of  such  a  process,  in  which  several  Ti3+  (Ce3+)  sites 
are  close  in  energy  [44,53].  Both  the  Cu— O  bond  lengths  (1.81  — 
1.91  A)  for  Ou — Ou  and  Ce-top  sites  are  close  to  the  experimental 
value  of  1.86  A  [54]  in  cubic  Q12O.  The  Cu— O  bond  length  (1.77  A) 
for  the  Ou  site  is  slightly  shorter  than  those  for  the  other  two 
adsorption  configurations  due  to  only  a  single  oxygen  atom  bonded 
with  the  Cu. 

Based  on  the  configurations  proposed  above,  we  calculated 
isosurfaces  of  excess  spin  density  for  a  Cu  atom  adsorbed  on 
CeC>2(100).  As  shown  in  Fig.  3,  the  isosurfaces,  which  are  denoted  by 
the  black  contour  lobes  in  each  of  Fig.  3(a)— (h),  represent  the  dif¬ 
ference  between  the  spin  up  and  spin  down  moments  at  given  Ce 
and  Cu  sites.  The  localized  spin  electrons  on  Ce  sites  in  Fig.  3  have 
/  ifxyz)  characteristics  [55],  The  contour  lobes  are  therefore 


indicative  of  charge  transfer  from  the  Cu  atom  to  one  of  the  Ce  sites 
(INN,  2NN  or  3NN)  listed  in  the  first  column  of  Table  1  due  to  the 
local  configuration  that  results  from  the  applied  distortion  and 
relaxation  (denoted  as  conf.).  Except  for  the  conf.  1  of  Ce-top  site 
where  two  2NN  Ce4H  ions  are  both  partially  reduced  (see  the 
smallest  jice  and  shortest  dce-o  values  in  Table  1),  all  other 
adsorption  configurations  involve  the  complete  reduction  of  NN 
Ce4~  ions,  as  indicated  by  magnetic  spin  moments  among  0.86  and 
0.96^b  which  are  close  to  the  ideal  moment  of  I^b  per  Ce3+  [56], 
Furthermore,  the  local  environments  of  the  Ce  ions  differ  from  one 
another  due  to  the  larger  size  of  the  Ce3+  ion  compared  with  the 
Ce4+  ion.  The  average  distances  between  surface  (or  subsurface) 
Ce3+  ions  and  the  neighboring  O  ions  for  the  respective  configu¬ 
rations  are  in  the  range  of  2.43—2.57  A  which  are  larger  than  the 
eight-fold  theoretical  and  experimental  bond  length  (2.38  and 
2.34  A)  of  Ce4+— O  in  cubic  CeC>2  (lattice  constants  of  5.49  and 
5.41  A),  but  are  comparable  to  the  average  experimental  Ce3+— O 
bond  length  of  2.455  and  2.542  A  in  cubic  (space  group:  IA-3,  six¬ 
fold,  206)  and  hexagonal  (space  group:  P3-M1,  seven-fold,  164) 
Ce2C>3  [57].  This  suggests  that  the  adsorption  is  accompanied  by  the 
reduction  of  Ce4+  to  Ce3+. 

For  the  most  energetically  favorable  configuration  of  Cu  at  the 
two-fold  Ou_Ou  site,  the  transferred  charge  localized  at  one  INN  Ce 
ion  leads  to  a  0.96jitB  spin  moment,  indicating  the  reduction  of 
Ce4+  — >  Ce3+.  The  corresponding  spin  density  isosurface  is  shown 
in  Fig.  3(d).  Fig.  3(e)  and  (f)  shows  isosurfaces  for  the  energetically 
less  favorable  conf.  2,  2NN  (the  electron  localized  on  the  Ce  ion  at 
the  subsurface)  and  conf.  3,  3NN  where  charge  is  localized  to  the/ 
orbitals  of  the  2NN  and  3NN  Ce  ions  relative  to  the  adsorbed  Cu. 
Fig.  3(a)— (c)  shows  spin  density  contours  at  Ou  sites  due  to  the 
associated  distortions  (dce-o(A))  listed  in  Table  1  as  well  as  the 
position  of  the  Cu  ion  after  optimization.  From  our  calculations,  the 
most  energetically  favorable  configuration  is  Cu  adsorption  at  the 
two-fold  Ou— Ou  site  with  the  electron  localized  at  the  first  nearest 
neighbor  (INN)  Ce  ion  (located  at  the  surface);  this  has  a  computed 
adsorption  energy  of  4.08  eV.  This  is  quite  different  from  the 
adsorption  behaviors  of  Cu  on  CeC>2(lll)  and  CeO2(110),  in  which 
the  adsorptions  at  the  three-fold  Od  site  with  the  excess  electron 
localized  at  the  third  nearest  neighbor  (3NN)  Ce  ion  and  at  the  two¬ 
fold  Ou— 0LI  site  with  the  electron  localized  at  the  second  nearest 
neighbor  (2NN)  Ce  ion  are  the  most  energetically  favorable  [15,17], 
On  the  other  hand,  for  the  same  adsorption  model,  the  adsorption 
of  a  Cu  atom  on  CeO2(100)  is  more  energetically  favorable  than 
those  on  Ce02(lll)  and  CeO2(110)  [15,17],  suggesting  that  the  more 
reactive  sites  exist  on  Ce02(100)  compared  with  those  on  Ce02(lll) 
and  CeO2(110).  Fig.  3(g)  shows  spin  density  contours  for  a  Cu  atom 
adsorbed  at  Ce-top  site  (conf.  1)  in  which  the  excess  electron  is 
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Table  3 

Adsorption  energies,  bond  lengths  and  magnetic  spin  moments  of  metal  atoms  for  the  adsorption  configurations  of  Au  atom  on  CeO2(100).  For  a  given  adsorption  site, 
configuration  numbers  (e.g.  conf.  2,  INN)  are  associated  with  the  distortion  applied  to  Ce— 0  distances  listed  as  dCe- o.  which  are  the  bond  lengths  between  the  Ce3+  and  the 
nearest  neighboring  0  ions.  The  undistorted  dCe- o  is  2.38  A.  £ads  is  the  computed  adsorption  energy,  dAu_0  is  the  distance  between  the  adsorbed  Au  and  the  nearest  0  ions  (an  0 
ion  for  the  Ou  site  and  two  0  ions  for  the  Ou-Ou  and  Ce-top  sites),  and  dAu-ce  are  distances  from  the  adsorbed  Au  to  INN,  2NN,  and  3NN  Ce  ions  (INN,  the  first  nearest  neighbor, 
2NN,  the  second  nearest  neighbor  and  3NN,  the  third  nearest  neighbor  Ce  ion  from  the  Au  adatom),  respectively,  following  geometry  optimization.  Note  that  integers  in  dAu-o 


(A)  and  dAu_Ce  (A)  denote  the  number  of  0  and  Ce  ions,  respectively,  within  the  distances  relative  to  the  adsorbed  Au  ion.  Magnetic  spin  moments  of  reduced  Ce  (in  bold)  and 
oxidized  Au  ions  are,  respectively,  nCe  and  nAu. 

Site  Eads  (eV) 

dAu-0  (A) 

dAu-Ce (A) 

PCeiPB) 

dce-O  (A) 

/“Au(/Ub) 

Ou  (conf.  1)  1.234 

1  X  1.970 

1  x  3.800 

0.496 

1  X  2.408 

0.000 

1  x  3.800 

0.495 

2  X  2.270 

4  x  5.375 

0.000 

1  x  2.592 

2  x  6.210 

0.000 

1  x  2.219 

2  x  6.679 

0.000 

1  x  2.422 

Ou  (conf.  2,  INN)  1.458 

1  x  1.970 

1  x  3.749 

0.000 

1  x  2.319 

0.000 

1  x  3.879 

0.956 

1  x  2.326 

1  x  5.315 

0.000 

1  x  2.337 

1  x  5.330 

0.000 

1  x  2.458 

1  x  5.447 

0.000 

1  x  2.519 

1  x  5.461 

0.000 

1  x  2.711 

1  x  6.231 

0.000 

1  x  6.236 

0.000 

1  x  6.488 

0.000 

1  x  6.888 

0.000 

Ou  (conf.  3,  2NN)  1.355 

1  x  1.973 

1  x  3.803 

0.000 

1  x  2.336 

0.000 

1  x  3.810 

0.000 

1  x  2.340 

1  x  5.276 

0.000 

1  x  2.410 

1  x  5.305 

0.000 

1  x  2.418 

1  x  5.327 

0.951 

1  x  2.599 

1  x  5.329 

0.000 

1  x  2.600 

1  x  6.162 

0.000 

1  x  6.172 

0.000 

1  x  6.627 

0.000 

1  x  6.655 

0.000 

Ou  (conf.  4,  3NN)  1.361 

1  x  1.974 

1  x  3.795 

0.000 

1  x  2.332 

0.000 

1  x  3.810 

0.000 

1  x  2.337 

1  x  5.365 

0.000 

1  x  2.385 

1  x  5.366 

0.000 

1  x  2.426 

2  x  5.373 

0.000 

1  x  2.545 

1  x  6.609 

0.000 

1  x  2.609 

2  x  6.228 

0.000 

1  x  6.642 

0.949 

Ou-Ou  (conf.  1 )  2.655 

2  x  2.027 

4  x  3.208 

4  x  0.277 

1  x  2.443 

0.000 

1  x  4.266 

0.000 

1  x  2.248 

2  x  5.762 

0.000 

1  x  2.274 

2  x  5.801 

0.000 

1  x  2.542 

4  x  6.317 

0.000 

1  x  2.189 

1  x  2.405 

Ou-Ou(conf.  2,  INN)  3.039 

1  x  2.021 

1  x  3.180 

0.000 

1  x  2.327 

0.000 

1  x  2.027 

1  x  3.184 

0.000 

1  x  2.346 

1  x  3.248 

0.000 

1  x  2.351 

1  x  3.261 

0.962 

1  x  2.487 

1  x  4.276 

0.000 

1  x  2.603 

1  x  5.766 

0.000 

1  x  2.661 

1  x  5.768 

0.000 

1  x  5.782 

0.000 

1  x  5.850 

0.000 

1  x  6.261 

0.000 

1  x  6.297 

0.000 

1  x  6.333 

0.000 

1  x  6.943 

0.000 

Ou-Ou  (conf.  3,  2NN)  2.846 

1  x  2.023 

1  x  3.216 

0.000 

1  x  2.399 

0.000 

1  x  2.024 

1  x  3.220 

0.000 

1  x  2.420 

2  x  3.230 

0.000 

1  x  2.426 

1  x  4.283 

0.853 

1  x  2.486 

1  x  5.796 

0.000 

1  x  2.487 

1  x  5.800 

0.000 

1  x  2.581 

1  x  5.826 

0.000 

1  x  2.572 

1  x  5.827 

0.000 

1  x  2.396 

1  x  6.326 

0.000 

1  x  6.328 

0.000 

1  x  6.333 

0.000 

Ou-Ou  (conf.  4.  3NN)  2.966 

2  x  2.025 

1  x  3.193 

0.000 

1  x  2.334 

0.000 

1  x  3.208 

0.000 

1  x  2.347 

1  x  3.217 

0.000 

1  x  2.388 

1  x  3.221 

0.000 

1  x  2.416 

1  x  4.280 

0.000 
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Table  3  (continued) 
Site 


Ce-top  (conf.  1 ) 


Ce-top  (conf.  2,  2NN) 


fads  (eV) 


0.944 


1.052 


dAu-o  (A) 


2  x  2.252 


1  x  2.192 
1  x  2.306 


4au— Ce  (A) 

1  x  5.760 
1  x  5.791 
1  x  5.800 
1  x  5.820 
1  x  6.261 
1  x  6.331 
1  x  6.352 
1  x  6.955 

1  x  2.795 

2  x  4.667 
2  x  4.669 


1  x  2.807 
1  x  4.581 
1  x  4.673 
1  x  4.671 

1  x  4.762 


MCc(Mb) 

0.000 
0.000 
0.000 
0.000 
0.940 
0.000 
0.000 
0.000 
0.000 
0.000 
2  x  0.523 


0.000 

0.000 

0.000 

0.000 

0.929 


dce-O  (A) 
1  X  2.558 
1  x  2.622 


1  x  2.211 

2  x  2.281 
1  x  2.433 
1  x  2.514 
1  x  2.576 
1  x  2.314 
1  x  2.315 
1  x  2.358 
1  x  2.543 
1  x  2.592 
1  x  2.622 


MAu(Mb) 


0.000 


0.000 


Distances  between  Au  and  reduced  Ce  are  labelled  in  bold. 


delocalized  at  the  2  s  nearest  neighboring  (2NN)  Ce  ions  and  the 
adsorption  energy  of  2.221  eV  is  the  lowest  among  all  the  consid¬ 
ered  adsorption  configurations.  Interestingly,  if  some  structural 
distortions  are  made,  the  conf.  1  at  the  Ce-top  site  will  evolve  into 
the  conf.  2  (see  Fig.  3(h)),  in  which  the  excess  electron  is  localized 
on  one  2NN  Ce4+  ion  and  the  adsorption  energy  increases  to 
2.638  eV. 

For  all  configurations  associated  with  each  adsorption  site,  the 
interaction  between  the  Cu  atom  and  CeO2(100)  surface  involves 
a  charge  transfer  from  the  adsorbate  to  the  substrate,  yielding  the 
oxidization  of  Cu  atom  and  the  reduction  of  a  Ce  ion.  Here  we 
take  the  most  energetically  favorable  adsorption  configuration 
(conf.  1,  INN  of  Ou — Ou  site)  for  an  example  and  analyze  the 
respective  partial  density  of  states  (PDOS)  of  Cu  atom,  pure 
CeO2(100)  and  Cu/CeC>2(100)  as  shown  in  Fig.  4(a)— (c),  respectively. 
In  Fig.  4(a),  there  are  occupied  spin-up  and  unoccupied  spin-down 
states  just  below  and  above  the  fF  which  correspond  to  the  Cu-4s 
state;  occupied  spin-up  and  spin-down  states  are  the  Cu-3d  state. 
For  pure  Ce02(100),  the  highest  occupied  valence  band  in  Fig.  4(b) 
results  from  0-2 p  and  some  contribution  from  the  Ce-4 f  suggest¬ 
ing  Ce-4//0-2p  hybridization  [46]  (herein  O  atom  refers  to  the  one 
which  bonds  with  Cu  atom  and  Ce  atom  refers  to  the  one  which  is 
reduced  upon  the  Cu  adsorption.).  We  can  clearly  see  that,  compared 
to  the  PDOS  of  the  isolated  Cu  atom  in  Fig.  4(a),  the  adsorption  of  Cu 
atom  at  an  Ou_Ou  site  (conf.  2  INN)  on  CeO2(100)  causes  the 
occurrence  of  a  new  state  just  below  £F  and  the  disappearance  of  Cu- 
4s  state  (see  Fig.  4(c)).  This  indicates  that  the  single  CU-4S1  electron 
has  been  transferred.  Simultaneously,  there  is  a  new  occupied  peak 
corresponding  to  the  Ce-4 f  states  for  the  Cu/Ce02(100)  system.  By 
comparing  the  PDOS  among  the  same  adsorption  configurations  of 
Cu  adsorbed  on  Ce02(100),  (110)  and  (111),  it  is  found  that  the 
Cu/Ce02(100)  exhibits  the  strongest  Cu-3d/0-2p  orbital  hybridiza¬ 
tion  interaction,  followed  by  Cu/CeO2(110)  and  Cu/Ce02( 111 ).  This  is 
also  consistent  with  the  relationship  among  the  adsorption  energies 
of  Cu  on  three  surfaces.  Moreover,  the  Cu-4s'  electron  localizes  on 
the  1 NN  Ce  ion,  resulting  in  reduction  of  the  Ce4+  ion  to  Ce3+  and 
oxidation  of  the  Cu  atom  to  Cu+.  This  observation  is  supported  by  the 
experimental  result  [58],  which  indicates  that  Cu+  species  are  al¬ 
ways  catalytic  centers  with  Ce02  playing  a  support  role. 

3.2.  Ag/CeO2(100) 

Following  the  same  methodology,  we  consider  the  same 
adsorption  sites  and  adsorption  configurations  as  in  the  case  of  Cu. 


Table  2  summarizes  the  adsorption  energies,  the  location  of  Ce3+, 
the  Ag— O,  Ag— Ce  and  Ce— O  bond  lengths,  and  magnetic  spin 
moments  of  Ag  and  Ce.  It  is  found  that  the  adsorption  at  the  two¬ 
fold  Ou— Ou  site  with  the  excess  electron  localized  at  the  first 
nearest  neighbor  (INN)  Ce  ion  (conf.  2,  INN)  is  the  most  ener¬ 
getically  favorable,  followed  by  those  at  the  Ce-top  and  Ou  sites.  The 
stability  sequence  for  the  Ag  adsorption  at  the  various  sites  on 
Ce02(100)  surface  is  similar  to  that  for  the  Cu  adsorption.  For  three 
adsorption  configurations  of  Ag  at  the  Ou— Ou  site,  the  computed 
fads  is  lower  by  1.642  (conf.  2,  INN),  1.653  (conf.  3,  2NN)  and  1.715 
(conf.  4,  3NN)  eV  than  those  of  Cu,  whereas  the  respective  ad¬ 
sorptions  of  Ag  at  the  Ce-top  and  Ou  sites  are  less  stable  by  0.968 
and  0.755—0.939  eV  than  those  of  Cu.  This  is  consistent  with 
relationship  between  the  Ag— O  and  Cu— O  bond  lengths  for  the 
same  adsorption  model.  Among  all  three  considered  sites,  the  Ag— 
O  bond  lengths  are  2.062—2.326  A  whereas  the  Cu— O  bond  lengths 
are  1.768—1.907  A  Note  that  for  the  adsorption  at  the  Ou— Ou  or  Ou 
site,  the  Ag— O  distances  of  2.062—2.078  A  are  in  close  accordance 
with  the  experimental  value  of  2.05  A  in  cubic  Ag20  [59], 

In  the  similar  situation  to  Cu,  the  adsorption  of  Ag  atom  at  dif¬ 
ferent  sites  has  different  configurations.  For  the  Ou  or  Ou— Ou  site, 
there  are  three  adsorption  configurations  in  which  the  Ag+— Ce3+ 
distances  are  different  (see  the  fourth  column  in  Table  2).  They 
are  separated  by  the  adsorption  energy  differences  of  0.058— 
0.232  eV  and  correspond  to  the  locations  of  Ce3+  at  the  INN,  2NN 
and  3NN  relative  to  the  Ag,  respectively,  in  which  the  INN  and  3NN 
are  in  the  surface  (the  second  atomic  layer)  whereas  the  2NN  is  in 
the  subsurface  (the  fourth  atomic  layer).  For  the  Ce-top  site,  the 
slight  structural  distortion  based  on  the  conf.  1,  in  which  the  excess 
electron  is  delocalized  at  the  2  s  nearest  neighboring  (2NN)  Ce  ions, 
results  in  the  increase  of  the  adsorption  energy  by  0.312  eV  and  the 
complete  localization  of  the  excess  electron  at  a  single  2NN  Ce  ion. 

According  to  Table  2,  we  can  see  that  upon  the  adsorption  of  Ag 
on  Ce02(100),  there  is  always  one  electron  transferred  for  all  the 
configurations.  In  most  cases,  the  transferred  electron  localizes  at 
one  Ce4 1  ion,  making  it  reduced  to  Ce3+,  companying  the  oxidation 
of  Ag  atom  into  Ag+  ion.  Fig.  5  compares  the  PDOS  of  a  Ag  atom 
(Fig.  5(a)),  pure  Ce02(100)  (Fig.  5(b))  and  the  most  energetically 
favorable  Ou— Ou  configuration  (conf.  2,  INN)  for  Ag/Ce02(100) 
(Fig.  5(c)).  As  is  the  case  for  a  Cu  adatom,  there  is  a  new  occupied 
peak  corresponding  to  the  Ce-4 f  state;  the  occupied  Ag-5s  state 
completely  disappears  (Fig.  5(c)).  This  implies  that  the  Ag  adatom  is 
oxidized  to  Ag+,  whereas  the  corresponding  Ce4+  is  reduced  to 
Ce3+.  This  is  consistent  with  calculated  results  of  Ag  adsorbed  on 
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(a)  Ou  (conf.  1 ) 


(b)  Ou  (conf.  2,  1 NN)  (c)  Ou  (conf.  3,  2NN) 


(d)  Ou  (conf.  4,  3NN)  (e)  Ou-Ou  (conf.  1 ) 


(f)  Ou-Ou  (conf.  2,  1  NN)  (g)  Ou-Ou  (conf.  3,  2NN)  (h)  Ou-Ou  (conf.  4,  3NN) 


(i)  Ce-top  (conf.  1)  (j)  Ce-top  (conf.  2,  2NN) 

Fig.  3.  Isosurfaces  of  excess  spin  density  for  Cu  adsorbed  (a)  Ou  (conf.  1 ),  (b)  Ou  (conf.  2, 1  NN),  (c)  Ou  (conf.  3, 2NN),  (d)  Ou  (conf.  4, 3NN),  (e)  Ou-Ou  (conf.  1 ),  (f)  Ou— Ou  (conf.  2,  INN), 
(g)  Ou — Ou  (conf.  3,  2NN),  (h)  Ou-Ou  (conf.  4,  3NN),  (i)  Ce-top  (conf.  1),  (j)  Ce-top  (conf.  2,  2NN)  on  CeO2(100).  Isosurfaces  correspond  to  the  location  of  Ce3+  ions  after  initial 
structural  distortions  followed  by  geometry  optimization.  Values  of  isosurfaces  are  0.002e  A-3  in  (a)— (j). 
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Fig.  4.  Partial  density  of  states  (PDOS)  of  (a)  Cu  atom,  (b)  pure  Ce02(100)  and  (c)  Cu 
adsorbed  on  Ou-Ou  site  (conf.  2,  INN).  Zero  energy  corresponds  to  the  Fermi  energy, 
£ F.  Upper  and  lower  panels  represent  spin  up  and  spin  down,  respectively. 


Fig.  5.  Partial  density  of  states  (PDOS)  of  (a)  Ag  atom,  (b)  pure  CeO2(100)  and  (c)  Ag 
adsorbed  on  Ou-Ou  site  (conf.  2,  INN).  Zero  energy  corresponds  to  the  Fermi  energy, 
£f.  Upper  and  lower  panels  represent  spin  up  and  spin  down,  respectively. 


Ce02(lll )  and  (110)  surfaces,  in  which  the  adsorbed  Ag  is  oxidized 
to  Ag+  [15,17].  Experimentally,  X-ray  photoelectron  spectroscopy 
(XPS)  analysis  of  Sarode  et  al.  [60]  indicates  that  the  high  catalytic 
activity  of  Ag  in  Ag/Ce02  is  attributed  to  the  presence  of  Ag+  species 
dispersed  on  CeCE  surface. 

3.3.  Au/Ce02(100) 

Au/CeC>2  catalyst  exhibits  good  catalytic  activities  in  several 
important  oxidation  reactions  [9,61—63]  and  Au/Ce02(lll)  and 
Au/CeO2(110)  have  been  extensively  investigated  [16,64—66]. 
However,  far  less  is  known  about  Au/CeO2(100).  Table  3  summa¬ 
rizes  the  adsorption  energies,  the  location  of  Ce3+,  the  relevant  Au— 
O,  Au— Ce  and  Ce— O  bond  lengths,  and  magnetic  spin  moments  of 
Ag  and  Ce  for  all  the  adsorption  configurations.  Similar  to  the  cases 
of  Cu  and  Ag,  the  adsorption  configuration  at  the  two-fold  Ou— Ou 
site,  in  which  one  INN  Ce4+  relative  to  the  Au  atom  is  reduced  to 
Ce3+  with  a  magnetic  moment  of  0.96^  and  the  Au  atom  is  oxi¬ 
dized  to  Au+,  is  the  energetically  most  favorable  (£adS  =  3.039  eV). 
Different  from  those  of  Cu  and  Ag  at  the  Ce-top  site,  the  second 
stable  configuration  of  Au  on  CeC>2(100)  is  at  the  Ou  site. 

Similar  to  the  cases  of  Cu  and  Ag,  there  are  different  configu¬ 
rations  of  Au  on  CeC>2(100)  with  the  excess  election  located  at 
various  neighboring  Ce  ion  and  they  are  separated  by  the  adsorp¬ 
tion  energy  differences  of 0.006—0.193  eV  (see  the  fourth  column  in 
Table  3).  Note  that  the  conf.  3,  2NN  and  conf.  4,  3NN  at  the  Ou  site 
have  the  same  adsorption  energies  (1.355  vs.  1.361  eV).  For  the 
configuration  at  the  Ce-top  site,  as  compared  with  the  delocalized 
solution,  an  0.108  eV  increase  of  the  adsorption  energy  for  Au, 


which  is  derived  from  the  slight  structural  distortion  and  the 
resultant  complete  localization  of  excess  electron,  is  smaller  than 
0.417  eV  for  Cu  and  0.312  eV  for  Ag. 

Fig.  6  shows  the  PDOS  of  an  isolated  Au  atom  (Fig.  6(a)),  pure 
Ce02(100)  (Fig.  6(b))  and  the  most  energetically  favorable  Ou— Ou 
(conf.  2,  INN)  configuration  (Fig.  6(c))  for  Au/CeO2(100).  Fig.  6(c) 
shows  that  the  adsorption  of  Au  atom  causes  the  occurrence  of 
a  new  occupied  peak  corresponding  to  localized  Ce-4 f  states.  The 
occupied  Au-6s  state  in  Fig.  6(a)  disappears  in  Fig.  6(c),  indicating 
that  Au  indeed  loses  an  electron.  Experimental  and  theoretical 
studies  [67—70]  proved  that  Au+  could  be  considered  as  an  active 
site  for  the  adsorption  of  CO  in  the  WGS  reactions.  Also,  experiment 
[71  ]  carried  out  by  Qian  et  al.  revealed  that  the  existence  of  Au+  in 
the  Au/Ce02  system  has  promoted  the  formation  of  surface  va¬ 
cancies  which  play  a  key  role  in  the  redox  reaction. 

3.4.  Discussion 

For  all  three  considered  adsorption  sites  on  Ce02(100),  Cu,  Ag 
and  Au  atoms  exhibit  the  similar  adsorption  characteristics.  The 
adsorption  at  the  two-fold  Ou— Ou  site  is  the  most  energetically 
favorable.  At  the  Ou— Ou  or  Ou  site,  there  are  three  different  con¬ 
figurations  with  the  excess  electron  (from  the  M  atom)  at  various 
neighboring  Ce  ions  (denoted  as  INN,  2NN  and  3NN).  Their  stability 
sequence  follows  £ads(l  NN)  >  £ads(3NN)  >  £a(js(2NN)  for  Cu  and  Au. 
However,  the  case  of  Ag  is  an  exception,  in  which  the  stability 
sequence  is  £a(js(3NN)  >  £a(js(lNN)  >  £ads(2NN)  at  the  Ou  site  and 
£ads(lNN)  >  £ads(2NN)  >  £a[|s(3NN)  at  the  Ou— Ou  site.  Interestingly, 
the  conf.  3,  2NN  and  conf.  4,  3NN  at  the  Ou— Ou  site  have  almost  the 
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Fig.  6.  Partial  density  of  states  (PDOS)  of  (a)  Au  atom,  (b)  pure  Ce02(100)  and  (c)  Au 
adsorbed  on  Ou-Ou  site  {conf.  2,  INN).  Zero  energy  corresponds  to  the  Fermi  energy, 
Ef.  Upper  and  lower  panels  represent  spin  up  and  spin  down,  respectively. 

same  adsorption  energies  although  the  2NN  and  3NN  Ce  ions  are 
located  in  the  subsurface  and  surface,  respectively.  For  the  Ce-top 
site,  there  are  two  configurations  in  which  the  excess  electron  is 
delocalized  at  the  two  2NN  Ce  ions  and  localized  at  a  single  2NN  Ce 
ion,  respectively,  and  the  latter  has  the  higher  adsorption  energy 
than  the  former. 

There  are  some  differences  among  the  adsorption  behaviors 
of  M  (M  =  Cu,  Ag,  Au)  on  Ce02(lll),  Ce02(110)  and  CeO2(100) 
that  require  clarification.  For  each  given  site  on  Ce02(lll)  and 
Ce02(110),  the  adsorption  configuration  with  the  excess  electron 
located  at  the  surface  Ce  ion  away  from  the  M  atom  is  energetically 
favorable.  The  opposite,  however,  is  true  for  the  M  atom  adsorption 
on  Ce02(100)  with  the  exception  of  the  adsorption  of  Ag  atom  at 
the  Ou  site.  In  addition,  it  is  worth  mentioning  that  the  adsorption 
of  Cu  atom  on  Ce02(110)  can  result  in  Cu2+  [17],  an  oxidation  state 
that  is  not  predicted  by  DFT  calculations  of  Cu  adsorption  on 
Ce02(lll)  and  Ce02(100).  The  orders  on  the  adsorption  energy  of  M 
(M  =  Cu,  Ag,  Au)  on  a  given  low-index  surface  follow  Cu  >  Au  >  Ag 
for  Ce02(100)  and  Ce02  (110)  and  Cu  >  Ag  >  Au  for  Ce02(lll).  On 
the  other  hand,  M/Ce 02(100)  (M  =  Cu,  Ag,  Au)  system  exhibits 
the  largest  adsorption  energy,  followed  by  M/Ce02(110)  and 
JVf/Ce02(lll),  which  is  consistent  with  the  relationship  among  the 
surface  energies  of  CeO2(100),  (110)  and  (111)  surfaces.  The  most 
reactive  surface,  (100),  [22—32,72]  provides  more  reactive  sites 
than  the  (111 )  and  (110)  surfaces  for  interaction  with  M{M  =  Cu,  Ag, 
Au).  There  are  compelling  experimental  results  that  support  this 
conclusion  [33,34]. 

The  importance  of  investigating  single  metal  atom  on  Ce02 
surface  arises  from  a  number  of  reasons,  including  (i)  nanocluster 


growth  proceeds  by  the  anchoring  of  one  atom  on  Ce02  surface, 
which  determines  future  growth  pathways  of  the  nanoclusters,  and 
(ii)  the  interaction  of  the  metal  atom  with  Ce02  surface  can  be 
probed  by  a  number  of  experimental  techniques  to  investigate  the 
fundamental  aspects  of  the  electronic  structure  of  interface  be¬ 
tween  metal  atom  and  Ce02  surface.  It  must  be  emphasized  that 
the  localization  of  the  excess  electron,  both  the  resultant  increase  of 
the  adsorption  energy  and  the  occurrence  of  multiple  adsorption 
configurations  with  the  excess  electron  localized  at  different 
neighboring  Ce  ion  are  attributed  to  the  effects  of  the  slight 
structural  distortions.  For  example,  for  the  undistorted  configura¬ 
tions  of  Cu,  Ag  and  Au  adsorbed  at  the  Ce-top  site,  the  excess 
electron  is  delocalized  at  two  2NN  Ce  ions  and  configurations  are 
not  stable  with  the  smaller  energies  of  0.108—0.417  eV  compared 
with  the  distorted  cases.  Similar  conclusion  has  been  reached  in  our 
previous  studies  of  M  adsorption  behavior  on  Ce02(lll)  [15]  and 
(110)  [17].  Also,  in  recent  studies  of  charge  transfer  behavior  caused 
by  the  removal  an  O  atom  in  bulk  Ce02  [42,73]  and  on  Ce02(lll) 
[44,45],  CeO2(110)  [74]  and  TiO2(110)  [53],  it  was  found  that  there 
exist  multiple  configurations  with  two  excess  electrons  localized  at 
different  Ce3+  ions.  As  mentioned  in  previous  sections,  Ce— O  dis¬ 
tances  in  each  configuration  detailed  in  Tables  1—3  were  slightly 
adjusted  to  lower  the  symmetry  of  the  corresponding  M/Ce O2(100) 
configuration  prior  to  geometry  optimization.  This  provided  the 
local  environment  that  was  conducive  to  the  reduction  of  Ce4+  to 
Ce3+  following  optimization.  Experimentally,  these  slight  structural 
distortions  can  result  from  the  introduction  of  impurities  or  by 
manipulation  of  the  crystal  morphology  during  catalytic  material 
preparation  [51,52].  The  combinations  of  experimental  result  with 
our  first-principles  modeling  deepen  the  understanding  of  the 
catalytic  properties  of  Ce02  surface  supported  metal  clusters  and 
could  also  be  expanded  to  other  important  metal/oxide  systems. 

4.  Conclusions 

Adsorption  mechanisms  of  M  (M  =  Cu,  Ag,  and  Au)  atom  on 
CeO2(100)  surface  have  been  investigated  using  the  first-principles 
density  functional  theory  calculations  considering  the  on-site 
Coulomb  interaction  within  the  PBE  +  U  scheme.  The  system¬ 
atical  analysis  of  the  adsorptions  energies,  atomic  and  electronic 
structures  reveals  that  the  order  on  the  adsorption  energy  follows 
the  sequence  of  Cu  >  Au  >  Ag.  There  are  multiple  adsorption 
configurations  at  a  given  site,  in  which  the  configuration  with  the 
excess  electron  localized  at  the  first  nearest  neighboring  Ce  ion  to 
the  JVf  is  the  most  energetically  favorable.  The  adsorption  behavior 
of  M  on  Ce02(100)  is  compared  with  those  on  Ce02(lll)  and 
CeO2(110). 
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